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1. Introduction and main result 

In theoretical fluid mechanics, the investigation of limit cases in which the thickness of 
the flow domain is small compared to its other lengthscale(s) is a classical subject. In most 
situations, the simplified equations describing these limit cases are derived formally from the 
original problem by expansion with respect to the small parameter describing the ratio of 
the lengthscales. Although lots of work have been devoted to the investigation of the limit 
equations (lubrication equations or various so-called Thin Film equations), the question of 
justifying the approximation by comparing the solutions of the original problem to those of 
the corresponding limit problem is less studied. This is true in particular when a moving 
boundary is an essential part of the original problem. 

In the case of 2D Hele-Shaw flow in a thin layer, the only rigorous limit result known 
to us has been proved by Giacomelli and Otto [9j. Their approach is based on variational 
methods and can even handle degenerate cases and complicated geometries. However, the 
existence of global smooth solutions of the Hele-Shaw problem under consideration has to 
be presupposed, and the obtained result on the closeness to some solution of the Thin Film 
equation is in a relatively weak sense and technically rather complicated. 

It is the aim of the present paper to provide a justification of the same limit equation 
using quite different, more standard methods. Starting from a strictly positive solution of 
the Thin Film equation, we prove the solvability of the corresponding moving boundary 
problems for large times. If the initial shape is smooth, approximations to arbitrary order 
and in arbitrarily strong norms are obtained. Moreover, our approach is straightforwardly 
generalizable to a multidimensional setting. However, it is restricted to the nondegenerate 
case of strictly positive film thickness and to a simple layer geometry. 

More precisely, we consider 2D Hele-Shaw flow in a periodic, thin liquid domain (i.e. a 
thread), symmetric about the x-axis, with surface tension as sole driving mechanism for the 
flow (see Fig. [T]). 
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Figure 1. The considered geometry and basic notation 



This problem (in the unsealed version, with surface tension coefficient normalized to 1) 
consists in finding a positive function h G C^([0, T], C^(§)), § = M/[0,27r), and, for each 
t G [0,T], a function u defined on 

n{h) := {{x,y) eSxR\\y\ < h{x)} 

such that 

-Au = in n{h), 

u = -K{h) on dn{h), (1.1) 
dth + Vu- i-h', = onr+(/i), 



where r+(/i) := {(x, h{x)) \ x G S}, and 

K{h) 



h" 



[I + (/l02)3/2 



is the curvature of r+(/i). (Here and in the sequel, for the sake of brevity we write h 
instead of h{t) or h{-,t) when no confusion seems likely. Moreover, we identify functions on 
§ with functions on 90 using the pull-back along x i— )• {x,zizh{x)).) Note that u represents 
the normalized pressure in the Hele-Shaw cell, and, in view of Darcy's law, the first and 
third equation are the incompressibility condition and the usual kinematic free boundary 
condition. Moreover, a uniqueness argument shows that u is symmetric with respect to y, 
i.e. u{x,y) = u{x,—y) for all (x,y) G 0(/i), meaning that Uy = on 8 x {0}. Thus, this 
setting corresponds to the case when the bottom of the Hele-Shaw cell, which we take as 
being § x {0}, is impermeable. We refrain from discussing further modelling aspects and 
refer instead to the extensive literature on the subject, see e.g. [4J, Ch. 1. Well-posedness 
results for (jl.ip (in slightly different geometric settings and various classes of functions) have 
been proved in e.g. [Tl[8l[13] and in [5, 6j for non-Newtonian fluids. 

To consider thin threads we introduce a scaling parameter e, < e <^ 1, and rescale by 

X = X, y = ey, h = eh, u{x,y) = u{x,ey). 
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Then u is defined on and {h,u) satisfies the e-dependent problem 

-e'^Uxx-Uyy = mQ{h), 

u = -eK{£,h) ondn{h), (1.2) 

{ dth + V^^^y)U-i-h',e-Y = on T+{h), 
where 

h" 

K(e, h) : 



(l + e2(/i')2)3/2 

Expanding u and Ke formally in power series in e we obtain 

u = eui + e^us + O(e^), 

where in particular 

iui)x = -h!" , {u^)y = hh"" on r+(^). 

So 

dth + e{hh"')' = 0{e'^), 

and after rescaling t = e~^t, suppressing tildes and neglecting higher order terms we obtain 
the well-known Thin Film equation 

dtho + {hoh'^y = on §. (1.3) 

Observe that (II. 3p is a fourth order parabolic equation for positive /iq which degenerates as 
ho approaches 0. For a review of the extensive literature on this and related equations we 
refer to [TT|. In the modelling context discussed here, (jl.Sp is used in [2l[3] to study the 
breakup behavior of Hele-Shaw threads. Note that (jl.Sp and its multidimensional analogue 

dth + dw{hVAh) =0, (1.4) 

also appear in models of ground water flow [llj . 

In view of the time rescaling, ho should be an approximation to 

■= [t ^ h{e~'^t)], 

where h solves (jl.2p . i.e. /i^ solves (with an appropriate v and omitting tildes) 

-s'^Vxx-Vyy = inQ{he), 
< V = —K{e,he) ondQ{hs), (1.5) 

^ dthe + Vv{-K,e-^)^ = onr+(/i,). 

(Of course, (II. 5p can be obtained immediately from (II. 1|] by choosing the "correct" scaling 
u = eu, t = e~^i at once. That scaling, however, is in itself rather a result of the above 
calculations.) 

Our rigorous justification of the Thin Film approximation for (jl.ip will therefore consist 
in showing that for any positive initial datum h* (from a suitable class of functions), the 
common initial condition 

/lo(0) = /le(0) = h* 

implies existence and uniqueness of solutions to (jl.Sp and (11. 5|) (for all sufficiently small e) 
on the same time interval, and 

he ^ ho as e I (1.6) 
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(in a suitable sense). This will be made precise in our main result Theorem 11.21 below. 
Observe that this, in particular, implies that the existence time in the original timescale, i.e. 
for solutions of (II. 2[) . goes to infinity as £ becomes small. 

We will make use of the following preliminary, nonuniform well-posedness result for solu- 
tions to (|1.5p . It is not optimal with respect to the demanded regularity but this is not our 
concern here. 

Theorem 1.1. Let s > 7 be an integer and h* G H^^^iS) a positive function. Then we have: 

(i) Existence and uniqueness: Problem (jl.Sp with initial condition /ie(0) = h* has a 
unique maximal solution 



K G C([0, T,), n cH[o, T,), 

for some = Ts{h*) G (0, oo]. 

(ii) Analyticity: We have 

[{x,t) ^ /is(x,t)]|sx(o,T,) G C^{^ X (0,T,)). 

(iii) Blowup: // < oo then 



Existence and uniqueness of the solution on some time interval [0, Q^ih*)] can be proved by 
parabolic energy estimates and Galerkin approximations as in [13j (for a different geometry). 
The arguments given there can also be used to prove (iii) by a standard continuation argu- 
ment. More precisely, for initial values h* that satisfy h* > a > and < M with any 
positive constants a and M, the existence time Qsi^*) has a positive lower bound depending 
only on these constants. This implies the blowup result (iii). For a proof of (ii) we refer to 
the framework given in ^ which is applicable here as well. Essentially, analyticity follows 
from the analytic character of all occurring nonlinearities together with the translational 
invariance of the problem. 

We are going to state the main result now. It sharpens (|1.6p as it also gives the asymptotics 
of /le to arbitrary order n E N. This is achieved by imposing strong smoothness demands on 
the initial value (or correspondingly, on Hq). To avoid additional technicalities, we have not 
strived for optimal regularity results. 

Theorem 1.2. Let s,n gN, s > 10 be given. There is an integer (3 = f3{s,n) £ N such that 
for any positive solution 



depending on /iq only such that for all e G (0, Eq) 

(i) problem (11. 5p with initial condition hf,{0) = /io(0) has precisely one solution 




/lo e C([0, T], //^(S)) n cH[o, T], i/^~^(S)) 
of (11. 3p there are Eq = eo{s,n,ho) > 0, C = C{s,n,ho) and functions 

hi,...hn-ieC{[0,T],H%S)) 



K e C7([o,r],i/^-i(S))nci([o,r],i/^-^(S)) 



(ii) 



he - {ho + ehi + ... + e"-i/i„_i) 



c([o,r],H=(§)) 
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As expected, the functions /ii, /12, . . . satisfy the (hnear parabohc) equations arising from 
formal expansion with respect to e. For details, see Lemma 14.21 below. 

Both this theorem and its proof are in strong analogy to [10] . where the parallel problem 
for Stokes flow in a thin layer has been discussed. In this case, the limit equation, also called 
a Thin Film equation, is 

dth + ^div{h^VAh) = 0. (1.7) 

The remainder of this paper is devoted to the proof of Theorem II. 2i For this purpose, 
we first transform (jl.5p to a fixed domain and rewrite the problem as a nonlinear, nonlocal 
operator equation for he ■ Some scaled estimates are gathered in Section [2j In Section [3] we 
discuss estimates for our nonlinear operators, while Section S] provides the necessary details 
on the series expansions that are used. Finally, the proof of Theorem 11.21 is completed in 
Section \5\ 

Technically, the main difficulty in comparison with the unsealed problem is the fact that 
the elliptic estimates for the (scaled) Laplacian and related operators degenerate as e becomes 
small. To handle this, weighted norms are introduced, and estimates in such norms have to 
be proved. In particular, the coercivity estimate for the transformed and scaled Dirichlet- 
Neumann operator given in Lemma 13.31 will be pivotal. On the other hand, the loss of 
regularity can be compensated by higher order expansions and interpolation. Moreover, as 
in [10], the ellipticity of the curvature operator is crucial. Therefore, the corresponding 
problems with gravity instead of surface tension appear intractable by the approach used 
here. 



2. Scaled trace inequalities and an extension operator 

In the remainder of this paper we let F-i- := S x {±1} denote the boundary components of 
our fixed reference domain O := § x (—1, 1). We write H^{Q), H^{T±) for the usual L^— based 
Sobolev spaces of order s G R, while by definition i?^(90) := H^{r^) x H^{T^). Functions 
/ G ^^^^(r-i-) may be represented by their Fourier series expansions 

/(x,±l) = E/(p)e'^'^ xeS, 

with f{p) the p—th. Fourier coefficient of /. Consequently, the norm of / may be defined by 
the relation 

11/11^ := \ Y.^l + Py\Kp)? 

Similarly, functions w £ H'^{il.) may be written in terms of their Fourier series 

w{x,y) = "^Wp{y)e'P'', {x,y) G n, 

and the norm of w is given, for s > 0, by the following expression: 

lkll^= (e(ii^^'|i')'+n'mkiio)' 
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where I := (—1, 1). Given e E (0, 1) and s > 1, we introduce the fohowing scaled norms on 

\H\t ■■= \Mti + \\d2w\\ti + elkll^ , w G H'{n), 

which are equivalent to the standard Sobolev norm || • ||^ but, for e — )• 0, degenerate to a 
weaker norm. The scaling is indicated by the first equation in (jl.Sp . and the scaled norms will 
enable us to take into account the different behaviour of the partial derivatives of the function 
V in system (jl.Sp with respect to e when e — )• 0. To do this we first introduce an appropriate 
extension operator for functions / G H'^{dQ) and reconsider some classical estimates in the 
weighted norms. 

Given e G (0, 1), we set for simplicity 

:= {di,„d2,e) ■■= {edi,d2). 

Lemma 2.1. (a) There exists a positive constant C such that for all e G (0, 1) and 
w G H^{Q) the following Poincare's inequalities are satisfied: 

\\w\\^<c\\vM\o Hac = o; (2-1) 

\\w\\^ <Ce-^\\Ve'w\\o if [ wda = 0. (2.2) 



(b) There exists a positive constant C such that the trace inequality 

\t+l/2 



^llf + ^/^ll«^llf^i/2<C^lkllm,. (2.3) 



is satisfied by all e G (0, 1), t G N, and w G H^+'^{Q). 
(c) Given £ G (0, 1), there exists linear extension operators E± : H^~^^^'^(T±) — > H^~^^(^}) 
such that (-E'±/)|r± = /; E-i-f are even with respect to y and E± satisfy the estimates 

\\E±f\\ti,e < Ct{\\ff,^+e'/'\\ff^^^^^), (2.4) 

\\d2E^f\\-i/2 < C(\\ffj;,, + e'/^ff,%), (2.5) 



1/2 ^ ^ yWJ ll_l/2 WJ lll/2^ 

t G N, / G H^^^^'^{T-i-), with constants independent of e. 

Proof. The proof of (a) is standard while that of (6) is similar to that of [10; Lemma 3.1]. To 
show (c), set for f{x, 1) = ^p/(p)e*^^ 

E+f{x,y) ■.= ^Wp{y)e^P^, Wp{y) := y2e^lplfe^-i)/(p). 
p 

Then llu^pUg < |/(p)| ^-nd for p ^ 

"1 „ /■! „ \?f^\\2 



ll^pllo' = \m\' I ^ lyV^^IPK^'-i) dy < 2|/(p)|2 ^ ye^eWy^-^) dy < 
Similarly, for /c G N, 

< Ckil + ie\p\f''-'mp)\' < C,il + N2^-i)|/(p)| 



\m: 

2e\p\ 
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Consequently, 



<C||/||^', 



\d2E+f\\ 



n 2 

t+i 



< 



<C(||/||[+' + .||/||f;,// 



< CsWfW 



r+ 2 

i+l/2 • 



This proves ([231 . To verify ([23]), note that 

\\d2E+f\\'i\ ' = Y,{l + \p\')-'/'\w'^{l) 



and 



K(l)P<C(l + e2|p|2)|7(p)|2. 
This imphes the result for E^, the construction for is analogous. 



□ 



Using an appropriate smooth cutoff function, one can construct a linear extension operator 
E : H^+^/^{dn) such that {Ef)\dn = f, Ef is even with respect to y if 

/(•,—!) = /(•,!), and E satisfies the estimates 



af7 ^^l/2||f||9C 



t+1/2 



\\d2Ef\\%, < C[\\f\\%, + e'/'\\f\\f%), 
t G N, / G //*+^/^(3ri), with constants independent of e. 



(2.6) 
(2.7) 



3. Uniform estimates for the scaled and transformed Dirichlet problem 

In this section we prove uniform estimates for the solution of the Dirichlet problem con- 
sisting of the first two equations of (jl.Sp . by using the scaled norms defined above. To this 
end we first transform the problem (jl.Sp to the strip $7 by using a diffeomorphism depending 
on the moving boundary /ig. 

Let M. he Q ov T±, a > dimA^/2, t < a. We will repeatedly and without explicit 
mentioning use the product estimate 

\\ziZ2\\f^ < C\\zi\\t'\\z2\\-^, zi G H\M), Z2 G H^M). 

For the remainder of the paper, let s and sq be such that s,so + 1/2 G N, sq > 7/2, 
s > 2sq + 3. (For example, sq = 7/2 and any s > 10 is possible, cf. Theorem 11.21 ) For given 
a, M > 0, define the open subset Ug := U{s, M, a) of i7*(S) by 

Us := \h G H^§) : ||/i|L < M and min/i > a). 
Moreover, define the (trivial) maps (j)± : T± — > S by (j)±{x, ±1) = x. 

To avoid losing regularity when transforming the problem onto the fixed reference manifold 
0, we modify [101 Lemma 4.1] to obtain the following result: 
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Lemma 3.1 (Extension of h). There exists a map 

[h^h]eC{H''{S),H''+^/'^{n)), cj>3/2, 
with the following properties: 



(i) h is even, h 



h o (j)^, and ^2/1 



0; 



(ii) IfheUs and /3 G (0,1), then ■= [{x,y) 1 — ^ (x, y/i(x, y))] G Diff2(0, f](/i)). 

In a first step we use the diffeomorphism to transform the scaled problem (jl.Sp into 
a nonlinear and nonlocal evolution equation on §, cf. (|3.3p . Therefore, we note that if 
/te : [0, Te) — )• is a solution of the scaled problem (jl.Sp . then setting w := —v o ^h^-, we find 
that the pair {h^,w) solves the problem 

'' -A{e,h)w = in r^, 

w = K{£,h)o(f)± on r±, (3-1) 

dth = B{e, h)w on S, 

where A : (0, 1) xUg ^ C{H^~^/'^{Q), H'^^'^ /"^{Vt)) is the linear operator given by 

^(e, h)w := DiDiW, 

with 

and Oj, i = 1, 2 given by 



Di := e(9i + aid2) 
di{yh) 



D2 := 02^2, 
1 



a2 := 



d2iyh)' ' d2{yh) 

Furthermore, we define the boundary operator B : (0,1) x Us C{H''~^/'^{n), H'^^iE,)) 
by the relation 

B{e, h)w{x) := {-h'di{w o $"!) + e-^d2{w o $"!)) (x, /i(2;)) 

= [—h'{diw + aid2w) + e~'^a2d2w) {x, 1), x G S. 
It is not difficult to see that B{£, h) may be also written as 

'a,; 



B{e, h)w = e 



a2 



-Diw 



(f>,^ on S, 



(3.2) 



where ai,e := eai and 02,^ := 02- 

Given / G f/''*~^(S) and (e, /i) G (0, 1) xUs, we denote throughout this paper by w{e, h){f} 
the solution w of the Dirichlet problem 

j -DiDiW = in n, 

1 w = / o 0^ on r±. 

With this notation, problem (j3.ip is equivalent to the abstract evolution equation 

dth = F{e,h) (3.3) 

where we set 

F{e,h){f} :=B{e,h)w{e,h){f}, (3.4) 
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and the nonlinear and nonlocal operator T : (0, 1) xlAs ^ H^^^{B) is given by the relation 

T{e, h) := F{e, h){K{e, h)}, (e, h) G (0, 1) x Ks- (3.5) 

It will become clear from the considerations that follow that w, F, and depend smoothly 
on their variables, i.e. 

w E C°°((0,1) xU,,C{H'-^{S),H'-^/^{n))), 

F gC^ {{0,1) X Us, C{H''\S),H'''' {§))), (3.6) 
C7°°((0,l) xUs,H'-''{§)). 

We start by estimating w and its derivatives, and finish the section by proving estimates 
for the function F. Some of the proofs rely on the following scaled version of the integration 
by parts formula 

f ^ f ^ f (^i e , f CLi e r 

/ — vDiwdx = — / — wDivdx+ / — —vw da — / — —vwda, (3.7) 
Jn ^2 Jn 0,2 yr+ «2 Jr_ «2 

which is true for all functions v^w £ H^{Q). 

In order to prove estimates for the solution operator w of (j3.2p . we begin by analysing the 
solution operator corresponding to the same problem when both equations in l\3.2\i have a 
nonzero right hand side. As a first result we have: 

Proposition 3.2. There exist constants Eq, C depending only on Us o,nd sq such that for 
integer t G [l,so + l/2], / G H^^^''^{E), fi G H*'^{Q.), i = 0,1,2, and {e,h) G (0,eo) x the 
solution of the Dirichlet problem 



-DiDiW = /o + di^efi in Q, 
w = f o (j)± on T±, 



(3.8) 



satisfies 



Additionally, 



w 



t<C[Y^ + + e^/2||/||,_i/2 . (3.9) 



,1=0 



\\d2w\\%2 ^^(11 ll/^llo + \\h\\%2 + ll/llo + ^'/'ll/lli/2j • (3.10) 

Proof. Step 1. We show ()3.9p for t = 1. We will consider the case / = first. 
Using relation (13. 7p . we proceed as in |10l Lemma 3.2] and find 

f 1 f 1 f 1 f 1 

I ■= — DiwDiwdx = — / — wDiDiwdx= / — wfodx+ / — wdiefidx 
Jn «2 Jn a2 Jn «2 Jn ^2 

—wfodx- I —di^swfidx- I di^e ( —] wfidx, 
n «2 Jna2 Jn \a2j 



10 



BOGDAN-VASILE MATIOC AND GEORG PROKERT 



where we used integration by parts to obtain the last equahty. So 

2 



|f2||„„l|f2 



«IIO 



W 



1,£- 



On the other hand, 



1=0 



I > c WM 



c{\\Vsw\\^)\ 



provided e G (0,eo) and eo is sufficiently small (with a constant c independent of e). Using 
Poincare's inequality (j2.ip . the estimate follows. 

If / 7^ 0, we let z := t<; — / E H^{Q), where / := Ef. Then z = on dQ and z solves in O 
the equation —DiDiZ = /o + di^gfi, where 

fo ■■= fo - ed2aiDj - d2a2D2l h = h + Dj, /a = /a + eaiDi7+ aa^z/. (3.11) 



Using ()2.6p and the result for homogeneous boundary data, we conclude that (|3.9p holds with 
t = 1. 

Step 2. We show (f3TI)]l . Define 



and observe 



B{£,h)w:=±e —^DiW on r± 

02 



d2W = ±e^h{l + e^h'^)-^{B{£, h)w + h'diw) on r±. 



(3.12) 
(3.13) 



We start with the case / = again. Then diw = and thus it is sufficient to estimate 
e'^\\B{e,h)w\\^^2- ^^^^ purpose, pick ip e H^/'^{d^) and define u G H'^{Vl) to be the 
solution of the Dirichlet problem 

f -DiDiU = in 

1 u = tp on d^. 

Then, by the transformed version of Green's second identity, 

1 



B{e, h)w4> da 



an 



0-2 



n 0.2 
u ' 

«2. 



-DiDiWU , dx 



h + edifi + d2f2 , 
u ax 



n 



02 





f U 


dx + — f2dcr 


1 02 . 


7r_ a2 



—f2da. 

r+ a2 



Consequently, applying the result of Step 1 to u, 



an 



B{e, h)wip da 



<c [\\u\\i,Y.\\f^\\o +\\f2\t'i 



an 

/2\\'y\\i/2 



< c 



l^^^/2 + Ell/^llo 
i=0 / 



ao. 

1/2- 
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This implies ()3.10p for / = 0. To treat the general case, define /, fi and z as in Step 1. Then, 
by the preliminary result, 

\\d2w\\%,<\\d2z\\%, + \\d2Ef\\%,, 
<ELoII/«II? + 11/211^^/2- 

Elo\\M^ <Elo\\M^ + C\\Ef\\^^,, 

\\f2\\%2 < \\f2\\%2 + C {e'\\f\\f% + Wd^EffJl/,) , 

and the result follows from ()2.6p and (|2.7p . 

Step 3. We prove p.9p by induction over t. The case t = 1 has been treated in Step 1. 
Assume now ()3.9p for an integer t £ [1, sq — 1/2]. Differentiating both equations of (|3.8p with 
respect to x we find that diw satisfies 

-DiDidiw = /o + di^Ji in il, 
diw = f o (f)± on F-i-, 

where 

/o = 9i/o - di2ai^eDiW - diai^ed2ai,ed2W, 
fi = difi + ediaid2W, 
h = dif2 + diai^sDiW + ai^£diai^ed2W. 
Using this and the induction assumption, we conclude that 

\\diw\\Z <c(y.^ Wfit + WfWt + ■ (3.14) 

In order to estimate ||922'U^||^_ii we use the first equation of (j3.8p and the explicit represen- 
tation 

^(e, h)w := e^diiw + 2e^aidi2W + {e^a\ + a^)d22W + {e^diai + e'^aid2ai + a2d2a2)d2W (3.15) 
to obtain 

o /o + di,sfi - e^diiw - 2£^aidi2W - {e^diai + e^aid2ai + a2d2a2)d2W 
O22W = — 2 ' 

and see that 

\\d22w\\tl ^ ^ II/^II^ + W^^'^te + ■ (3.16) 

Combining (|3.14p . (j3.16p . the induction assumptions, and the relation 

Ikllm,. < C {\\w\\^^, + \\dM\t + \\d22w\\ti) 

yields the desired estimate for This completes the proof. □ 

Using this result we can additionally show that then 

\\d2w{e,h){f}\\l_^/2 < Ce^\\f\l^+,/2. (3.17) 

(Note that this involves a higher norm of /, but the constant involved in the estimate is of 
order e^.) 
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To show this, let G H^~'^{Q) be the extension of / given by cl){x,y) = f{x) and define 
z := w{e,h){f} - (p. Then 

{-DiDiZ = e'^diicj) in Q, 
z = on dQ, 

and by the unsealed trace inequality and ()3.9p with t = sq + 1/2 we get 

\\d2w\\l_y, = \\d,z\\l_,/, < + < Cs^\\dnm^^/2 

and therefore (j3.17p . In particular, this implies by the unsealed trace estimate 

\\d2w{e,h){fW]^_, < Ce'\\f\U,^y2. (3.18) 

Next, we prove a coercivity estimate for the scaled Dirichlet-Neumann operator F{e,h), 
{e,h) G (0.en)xUe. which will be a key point in the proof of Theorem [L2l Given </? G iJ^/^(§) 
and e > 0, we set 

IIV'lliAe := llv^llo + e"^''^||<^||i/2- 

Lemma 3.3. There exists a positive constant c such that for all (e,/i) G (0, eg) ^ CLnd 
^p £ H^/'^(§) which satisfy 



Is 



ipdx = 

JS 

we have 

{F{s,h){ip}\ip)L2^S)>cM\l/2,e- (3-19) 

Proof. Let w := w{e, h){ip} G H'^{Q), recall the definition of B{e, h)w from (I3.12p and observe 
that due to symmetry 

B{e,h)w{x,l) =B{e,h)w{x,-l) x G S. 

Using (|3.7p . we have 

{F{e,h){ip\\Lp) i2i^\ = I wB(e,h)w da = e^"^ / -^wDiwda 

e^^ f aie f ai ^ f 1 

' wDjW da / — —wDjW da = / — DjwDjW dx 



2 Jr+ ^2 2 Jy_ a2 2 Jq a2 

>ce-2 / iV.wl'^ dx = ce-^ {WVewWo)^ , 
Jo, 

cf. Proposition [321 From Poincare's inequality (12.2p together with (I2.3P we obtain the desired 
estimate. □ 



Now we prove estimates for the Frechet derivatives of the solution w = w{e, h){f} of (13. 2p 
with respect to h. The results established in Proposition 13.21 will be used as basis for an 
induction argument. 
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Proposition 3.4. Given k £ N, hi, . . . ,hk £ H^{§), and f G H^~'^{S), the Frechet derivative 
y^(k) := . . . ^ satisfies 

\\w^'''^\\t<C\\hiU--A\hk\\sMt-i/2. (3.20) 
for all integer t G [1, sq — 1/2]. Additionally, 

\\d2W^%%2<C\\h4so--A\hk\\sMi/2. (3.21) 
The constant C depends only on k, sq, andUg- 

Proof. We prove both estimates by induction over k. For /c = they hold due to Proposition 
13.21 Assume now (j3.20p . (j3.2ip for all Frechet derivatives up to order k. Differentiating (j3.2p 
{k + 1)— times with respect to h, yields that is the solution of 

' -D,DiW^^+^) = E EQa[+M(e,/i)K(i),.../i.(;+i)M'=-^) inn. 



o-GSi,_|_i l=Q 

w^'^+i) = on dn, 



(3.22) 



where w^^ = w^^ 'H^' ^)[^o-(«+2)) . . . ) ^(T{fc+i)]{/} aiid Sk+i is the set of permutations of 
{1,...,A: + 1}. 

We are going to define functions Fi in n such that the right hand side in (I3.22p i can be 
written as 

k 

The functions Fi are sums of terms to be specified below. For this purpose, we recall (|3.15p 
and consider the Frechet derivatives of the ocurring terms separately. 

(i) When differentiating e'^d22W we do not obtain any term on the right hand side of the 
first equation of (I3.22p . 

(ii) The terms on the right hand side of the first equation of (I3.22p which are obtained 
by differentiating 2e^aidi2W may be written as follows: 

e'a^l+'^di2W^'-^^ = di,e [eaf^'^d2W^''^^] - {e^af^'^) d2W^'-^ 

where af^^^ := af^^\h)[hf^(^i-^ , . . . /io-(i+i)]. The last term belongs to Fq, while the one 
in the square brackets belongs to Fi. 

(iii) When differentiating (e^af + a|)922W^ we obtain terms of the form 

{e^aj + a^)('+i)a22u;('-') = ^2,. [(e^a? + a',f+'^ d2w(''-'^] - d2ie'al + d2w'^''-'\ 

where (e^af + 02)*^'+^-' := 9^"^^(e^af + a2)(/i)[/io-(i)> . . . ^cr(«+i)]- The last term belongs 
to Fq while the expression in the square brackets belongs to F2. 

(iv) All terms corresponding to {e'^diai + e'^aid2a2 + a2d2a2)d2W are absorbed by Fq. 

Summarizing, we get 

r -AA?x'(^'+') = Fo + ^^,,Fi inn, 

\ = on an. 
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where 

k 

aeSk+i 1=0 

and 

with smooth functions /3/j and |7/o,j| G {0, 1, 2}, |7n,j|) l7/2,j| £ {0, !}• Fixing a and /, writing 
(^H ■= a/i[^o-(i)7 • • • 5 f^a{i+i)] ^-iid using the induction assumption we estimate 

WliOiW'- < 0||a/i||^^_3/2l|C2U'^ ^llt-l < '-'l|/i<7{l)llso+l/2 • • • ll'i<7{«+l)llso+l/2lF '\\t,e 

< C\\hi\\so . . . ||/lfe+l||sol|/||t-l/2- 

Thus 

\\Fi\\ti<C\\hi\U,...\\hk+^\\so\\f\\t-i/^ i = 0,1,2, (3.24) 
and (j3.20p (with k replaced by A; + 1) follows from p.9p . 
Similarly, 

\\ai202W'- ^||_i/2 ^ '-^ll"'2|lso-ll|0^2^^^ II -1/2 ^ '-^ II "«2 IUo_l/2 II 0*2 11-1/2 
< 11/^^(1)11^0+1/2 • • • ll'^a(/+l)IUo+l/2ll'^2W^ '\\-l/2 



< C'll^lIU • • • ll^fc+llUlU llt-1/2 

Therefore 

||F2||%<C||/li||,„...||/l,+i||,J|/||i/2, 

and (j3^ (with k replaced by A; + 1) follows from with t = 1 and (iXTO]) . □ 

We prove now an estimate similar to p.2ip which is optimal with respect to one of the 
"variations" /i^ (say hi). The price to pay here is a stronger norm for /. 

Proposition 3.5. Under the assumptions o f Proposition [3l4l we additionally have 

||a2^«('=)||-l/2 < C||/ii||3/2||/i2|Uo . . . ||/ifc|Uoll/l|so-l. 

The constant C depends only on k, sq, and Us- 
Proof. We show the more general estimate 

||u;W||?,, + ||a2«;('=)||-i/2 < C||/ii||3/2||/i2|Uo . . . \\hk\Uf\\so-i (3-25) 

by induction over k. For k = the statement is contained in Proposition 13.21 Assume now 
()3.25p for all derivatives up to some order k. We proceed as in the proof of Proposition 13. 4|, 
reconsider problem ()3.23p and have to show now 

\\Fi\\^, \\F2\\%, < C||/li||3/2||/l2|Uo • • • \\hk\U\f\\so-l. 

For this purpose, we fix a and / and estimate 

\\aiid2W^'-'^\^, ||a/292«;('-')||^^^/2. 

We have to distinguish two cases, depending on whether the argument hi occurs in the first 
or in the second factor. 
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Case 1: a^^ (1) < I + 1. Using Proposition 13.41 with t = sq — 1/2 we estimate 

^'lso-l/2,£ 



so • • • ll^fcllsoll/IUo-l- 
The product is taken over j G {1, ...,/ + 1} \ {cr^^(l)}. Similarly, 

\\ai2d2W^'-^^\\%, < C\\ai2\\!%\\d2w('''^C_, < C||az2||f ll^st.^^-') ||?„_3/2 

<C\\h,\\^ n ll^'^a)llS,+l/2 lk^'"'^llS,-l/2,. < C^I|/^lll3/2l|/^2|Uo • • • \\hk\U\f\\so-l- 

Case 2: cr^^(l) < / + 1. We apply the induction assumption and estimate 

<C||/ll||3/2l|/i2||so...||/ifc||.oll/IUo-l- 

Similarly, 

\\ai2d2W^'-'^\\%, < C\\ai2fs^,,\\d2W^''^^\\%, < C\M\i_,/2\\d2W^''Y^y2 
<C\\Kii)\\l+i/2-- -11^(1+^^^^^^^ 

The proof is completed now by carrying out the summations over a and / and applying 
Proposition EJ] to IK23\i . □ 

We recall (cf. 



He, h){f} = + h''){d2w{e, h){f})\r_, o - h'f. 

Applying the product rule of differentiation and product estimates as above we find from this 
and Propositions 13.41 and 13.51 

||F(-)(e, h)[h,,... M{/}|Ui/2 < Ce-^h,\U, . . . \\hk\\so\\f\\i/2, (3.26) 
||f("^)(6, h)[hi,... U{/}ll-i/2 < Ce-^h,h/2 . . . \\hk\\so\\f\\so-i (3.27) 

for all m e N, (e, /i) G (0, Cq) x ^s, / G -f/''*"^(S) , and hi, . . . hm € H^iS). Additionally, using 
(IXTHD . 

||F(e,/i){/}|U„_i<C||/||,„+3/2. 

In particular, we have 

||^(e,/i)||,o_i <C||/i||,,+7/2 <C, {e,h) e{0,eo)xUs. (3.28) 

The constants depend only upon Ug, sq, and m. 
Moreover, we obtain: 

Lemma 3.6. Given hi G //"*(§) and f G //*~^(S), we have 

\\F'{e,h)[hi]{f}\\_y2 < C||/ii||3/2||/||.o+3/2. (3-29) 
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Proof. For brevity we write w' := w' {e, h){f}. Differentiating ()3.2p with respect to h yields 

J -DiDiw' = fo + di^efi in ^, 
1 ui' = on d^}, 

where 

:=e'^dhai[hi]di2W + 2~^dh{e^al + al)[hi]d22'w, 
■.=£dhai[hi]d2W, 
■.=2-^dhie^al + al)[hi]d2W, 
■.=w{e,h){f}. 

By dSJ]) and (f3l7|) we have 

2 

\\d2w'\\y2 <\\w'\\Z <CY, WfiWl < C\\h,\\y2\\d2w\\^ < Ce^\\hih/2\\fK+3/2- 



k 
fi 
h 

w 



i=0 



The result follows easily from this. 



□ 



Next we give an estimate for the remainder term that occurs when curvature differences 
are linearized. 

Lemma 3.7. Let e G (0, 1) and h^heUsD i7*+3/2(S). Then 

dtr^ (K(e, h) - K{e,h)) - h)[{h - 7^)(^-i)] < C(l + ||7^||,+3/2)||/i - h\\s+i/2- (3.30) 

1/1 

The constant C depends only on Ug ■ 
Proof. By the chain rule, 

s-l 

dtr^nie, h) = K'{e, h)[h'^'-^^\ + ^ Cp,...p,K^^\e, h)[h^P^\ . . . , h^P'^], 

1=2 

with 1 < pi < . . . < pi, and pi + . . . + pi = s — 1. This also holds if we replace h by h. We 
subtract these identities and obtain 



{K{e, h) - k{£, h)) - k'{£, h) [{h - h) 



s-l 



1=2 



j=i / 
The terms on the right are estimated separately. One straightforwardly gets 

\\K^'\e, h)[hl,..., hl]h/2 < CWhlh . . . ||/ii_i||3||/lz||5/2 



(3.31) 



(3.32) 
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for all / G N, / > 1, hi E H^/^{S) and hi,... hi_i G H'^{S). Since Ug is convex, we additionally 
have 

'K^^){e,h)-K^'\e,h)) [hi,..., hi] 



1/2 



< 



n(i+^)(e, rh+{l- r)h)[h - h, hi, . . . , hi] dr < C\\h - hh\\hih ■ ■ ■ II Villall^; 



1/2 



l5/2- 



Applying these estimates to all terms in (j3.3ip and adding them up yields the result. □ 

Finally, we give a parallel estimate concerning the complete operator T. Using the invari- 
ance of our problem with respect to horizontal translations we obtain, as in [10], Eq. (6.8), 
the "chain rule" 

dtr^F{e, h) = F{e, h){dtr^K{e, h)} + J] C,,_p,^,F^^\e, h)[h^P-\. . . , h)}, 

k>l 

(3.33) 

h £Us sufficiently smooth. The sum is taken over all {k + l)-tuples {pi, . . . ,Pk+i) satisfying 
pi + ... +Pk+i = s-l and pi,... ,Pk > 1. 

Lemma 3.8. Additionally to Lemma \3.7\ assume £ £ (0,eo)- Define 

Ps{£,h,h) :=ar^ {F{e,h) - F{e,h)) - F{e,h){K' {e,h)[{h -h)'^'-^^]}. 

Then 

\\Psie,h,h)\\_y^ < Ce-\l + m,+3/2)\\h -h\\s+i/2- (3.34) 
The constant C depends only on Us ■ 



Proof. Observe that by density and continuity arguments, it is sufficient to show (j3.34p under 
the additional assumption that h and h are smooth. We infer from ()3.33p that Ps{£,h,h) = 
E" + + G, with 

E" :=F(e, h){d'-^K{£, h)} - F{£,h){d'^-\{£,h)} - F{e, h){K'{e, h)[{h - h^'^^ 

E^ ■.=F'{e,h)[h^'-^'^]{K{e,h)] - F'(e, 7^)[7^^'"^^]{K(e,7^)} 

and G := ^ Cvi,---,Pk+i-^pi,...,pk+v "^^^^^^ taken over all tuples satisfying additionally 

1 < Vk+i < s-2, and E^^^ =: E'' is given by 

E'^ :=F(^')(e, h)[h^P'\ /i(P^)]{^'=+iK(e, h)} - F^''\e,h)[h^^'\ . . . ,h^^'^]{dS:''+' K{e,h)}. 
We estimate E"" first and write E"" = Ef + E2, where 

Ef ■.=F{e, h){dtr^{K{e, h) - ti{e,h)) - K'{e, h)[{h - h)^'-^^]} 

E^ := {F{E,h) - F{e,h)) {dt^-^ie ,h)} . 
Invoking ()3.26p (with m = 0) and Lemma 13.71 we get that 

||£;?||_i/2 <Ce-'\\drH^ie,h)->^{e,h))-K'{e,h)[{h-h)^'-'^]\\i/2 

<Ce-\l + \\h\Us/2)\\h-h\Ui/2. 



(3.35) 
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In order to estimate E2, we write 

El= [ F'{e,rh + {l-r)h)[h-h]{d'^-^K{e,h)}dr, 
Jo 

and using (j3.26p . with m = 1, yields 

\\En-i/2<e-^\\h-MsM\s+3/2<Ce-^\\h\U3/2\\h-h\\s. (3.36) 
Similarly, we decompose = + E2 + E^, where 

E', := {F'{e,h) - F'{e,h)) [h^'~\^{e,h)} 

= / F"{e,rh + {l-r)h)[h-h,h''''\n{e,h)}dr, 
Jo 

E\ :=F'(e,/i)[7I^'"^^]{K(e,/i) -^(e,/^)}, 

E\ :=F'(e,/i)[/i(^-i) - li'"'^\{K{e,K)\. 
The estimate ()3.27p with m = 1 and m = 2, respectively, yields 

\\E\\Ui2<Ce-^\\h-^\,,\^^'-'-\\^l^^^^^^ 

< Ce-'\\h^'-'\\y2Me,h) - K{e,h)\\so < Ce-'mU+y^Wh -h\\s, 

\\El\U/2 < Ce-V'^ -h^'~'''h/2Me,h)\U, < Ce~^h -h\U,/2. 
To estimate G, we proceed similarly and decompose E"^ = Ei + E2 + E^, with 

El := h) - F^'HeJ)) [h^''\ . . . , 7^(^'=)]{aS'=+V(e,7I)} 

= / F^''+^\e,rh + {l-r)h)[h-h,h^'''\...,h'^''^]{d^''+'K{e,h)}dr, 
Jo 

k 

1=1 

We distinguish two cases. 

Case 1. Suppose first that Pk+i ^ Pj for all 1 < j <k. Then 

s-l 

Pk+i<s-2, pi, . . . < < s - So - 2, 

by the choice of s. Choosing m = A; + 1, we infer from Pk+i > 1 that k + 1 < s — 1, and 
together with relation ()3.26p we find 

\\Et\\_y2 < Ce-'\\h-h\Uh^''^\\so . . . ||7^^"=^||solHI.+i/2 < Ce-'\\h\U,/2\\h-h\\s, 
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while, for m = k, the same relation implies 

11^211-1/2 < Ce~'^\\h - /i||s+i/2, 

\\El\\_y2<Ce~^h\U,/2\\h-h\\s- 

Case 2. Due to symmetry, we only have to consider the case whenpi > pj, for all 1 < j < k+1. 
Then 

s - 1 

Pi < s-2, P2,--- ,Pk+i < < s - so - 2, 

and (|3.27|) with m = k + 1 and m = k, respectively, yields 

\\Ef\\^i/2<Ce-^h-h\\s. 
This completes the proof. □ 

4. Approximation by power series in e 

In this section we construct operators J-k and functions t i— )• hgr^kit) such that, in a sense 
to be made precise below, 

T{e,h) = +0(e^+i), 

and h^^k is an approximate solution to (j3.3p . Formally, the construction is by expansion with 
respect to e near 0, i.e., Fk{£-,h) and /le^jt are polynomials of order k in e. In lowest order 
A: = 0, we will recover the Thin Film equation (jl.3p . As this construction involves a loss of 
regularity that increases with k, we will have to assume higher smoothness of h. 

Fix /c G N and let si > s + A; + 15/2, with s as before. In this section, we will assume 
h G and all constants in our estimates will be independent of h. 

We start with a series expansion for iL'(e, 

Lemma 4.1. For p = 0, 1, . . . , A; + 2 there are operators 



w 



such that 

n 

k+2 



w{e,h){f}-Y,ePw^''Hh){f} 

p=0 



< Ce''+^fU-2. 

s+3/2 



In particular, 

-1 



= fix), w^''\h){f}{x,y) = f"{x) / -:^7—;dT, ^x;^ = ^ for p odd. 



ly a2(x,T) 

Proof. Recalhng (j3.15p we have A{e, h) = 5o(/i) + e'^S2{h) with 
So{h) -.=0^822 + 0202,2^^2, 

S2{h) ■■=dfi + 201(9^2 + ald22 + (ai,i + ai«i,2) ^2. 
The terms ?yfPl(/i){/} are determined successively from inserting the ansatz 



(4.1) 



k+2 

w{e,h){f} = Y,ePw^PHh){f} + R 

p=0 
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into 



iSo{h) + e^S2{h))w{e,h){f} = in Q, 

wie,h){f} = /o(/>± on r±, 
and equating terms with equal powers of e. Thus we obtain 



5o(/i)«;[' 



in n, 

f o(t)± on r± 



5o(/i)it;W 



w 



[1] 



in Q., 
on da. 



and further 



J 5o(/i)i«[P+2] = -52(/i)u;[P] in 0, 

\ ^0 on 517, 

p = 0, . . . , /c. Observe that the general problem 

{SQ{h)u = G in 0, 

ti = g o (t>± on F-i- 

(with (7 and G even) is solved by 



u{x,y) = g{x) 
and for this solution we have 



1 



a2{x,TY Jq 



G{x, s) ds dr, 



u 



<C{\\g\\t + \\G\\^), 



t £ [s + 3/2, si — 2]. All statements concerning the mapping properties and the explicit form 
of the w^P'^ follow from this. To estimate the remainder, observe that 

J A{e,h)R = -e'=+352(/i)u'[^+il(/i){/} -e'=+%(/i)u;['=+2](/i){/} in n, 

\ R = on dQ. 

The estimate follows from Proposition 13.21 with s replaced by si and t = s + 5/2. 

Recall, furthermore, that 

B{e,h) = e-^B^^^ih) + B^^^h) 

where 



□ 



are given by 

B^'^\h)w = h-\d2w)\r+ o 0-1, B^^\h)w = -h\diw)\r+ o 0;^ + h-\h')\d2w)\r+ o 
By Taylor expansion around e = it is straightforward to see that there are functions 
^[p] ec~(ZY,,,/7^i-2(S)), p = 0,...,k + 2, 

such that 

k+2 



p=0 



< Ge 



fc+3 



si-2 



In particular. 



1 (h) = h" and = for p odd. 
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In view of (13. 4|) , (|3.5p we define 

k+2 

p=0 j+m+l=p 

j E {0, 2}. As all terms corresponding to p = and p = 1 vanish, this is indeed a polynomial 
in £ and 

J-fc G C°°([0,1) xU.^H^S)). 

In particular, 

To{£,h) =Tk{0,h) = 13^^^{h)w^'^\h){K^^^{h)} + B^'^\h)w^^^{h){K^^\h)} = -{hh'")' 

(cf. ([13])). 

It is straightforward now to obtain 

\\F{e,h)-Ft,{eM\s<CE^+\ (4.2) 

To construct the approximation h^^k we start with an arbitrary, sufficiently smooth, strictly 
positive solution /iq of the Thin Film equation (II. 3p and successively add higher order cor- 
rections. We closely follow [101 Lemma 5.3] here. Fix T > 0, h* £ Us, and set for brevity 

T ■=k + 15/2, 

S2 := S2{k, s) -.= 3 + [k/2]{T - 4) + r + 1, 

Va ■■= {H e C{[0,T],Usr]H''{§)) nC\[0,T],H''-\E)) \H{0) = h*}, a>s. 
Let ho £ be a solution to ()1.3p . Observe that 

-FfceC°°([0,l) x^/,,//'^-^(§)), a€[s + T,S2]. (4.3) 
Furthermore, for t G [0, T], the linear fourth order differential operator 

A := A{t) := 8^0, ho{t)) = [h^ {hh'^{t) + ho{t)h"'y] 
is elliptic, uniformly in x and t. 

Lemma 4.2. Fix /iq as above. There are positive constants CLi^d C and functions 
he^k G Vs+4, e G [0,eo), that satisfy 

he,k{t)dx= I ho{0)dx and \\dthe,kit) - F{e,he,km\s < C£''+\ te[0,T]. (4.4) 

Proof. We construct /i^^^ by the ansatz 

he,k ■■= ho + ehi + . . . + e^hk, 

where for p = 1, . . . ,k, hp is recursively determined from ho, ... , hp-i as solution of the fourth 
order linear parabolic Cauchy problem 

dthp = ^^(-^fc(e, ^£,A:)|£=o = ^hp + Rp, 
hp{0) = 0, 

where Rp = [t -Rp(0] is a finite sum of terms of the form 

m 

did'^Tk{o,ho)[h,„...,h,j, i<ji<p-i, i + Y,j,=p. 

i=l 
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(At this point, the expression -^{Tki^, ^e,fe)|e=o should be understood in the sense of formal 
expansions only. It will be justified below.) Note that ^'^^^(O, /iq) = for / odd, and therefore 
also hp = for p odd. Fix 9 £ (0, 1/4). We will show by induction that 

hp G C^([0,T],F--(S)) n C'+\[0,T],H"''-\§)), ap:=S2-l- |(r - 4). 

For p = 0, this follow from our assumptions by a standard interpolation argument. Suppose 
now this is true up to some even p < k — 2. By (j4.3p we find 

Rp+2 G C^([0,T],F'^--^(S)), Rp{0) G H^'-^S) 

and by standard results on linear parabolic equations (cf. e.g. [I2l Prop. 6. 1.3]) 

hp+2 G C\[0,T], H^P'^+'^iS)) n 0^+^^ {[0,T], H^p-^ {§)). 

Therefore, by our choice of S2, 

h,^k G c([o,r], //«+-(§)) ncH[o,r],/7^+--4(§)). 

If e G [0,eo)) ^0 sufficiently small, this implies h^^k{t) G Us+r-, and thus, by Taylor's theorem 
applied at e = to e i— ?■ dth^ ^ — ^k{^-, ^e,fc); 

\\dtK^k-me,he,k)\\s<C£''+^ 

Consequently, we get (14. 4p from this and (14. 2p . 

Finally, for all h € IAs+t we have J^J-{e, h) dx = 0, cf. [5^ Lemma 3.1] and \1?>\ Lemma 1]. 
Therefore, by (j4.4p . jg dth^^^ dx = 0{e^^^). This implies f^hp dx = 0, p = I, . . . , k, and thus 
the lemma is proved completely. □ 

5. Proof of the main result 
Let T, Us, h* as in the previous section and fix T' £ (0,T]. Let 

h, e c{[o,r],Us) nC^{[o,T'],H'-^{8)) 

be a solution of (|3.3p with hs{0) = h*. For given, sufficiently smooth ho solving (|1.3p . we 
denote by h^^k the function constructed in Lemma 14.21 The following energy estimates are 
the core of our result. 

Proposition 5.1. 

(i) Fixk £ N and a solution /iq G Vs^(^k) of ()1.3p . There are constants C and depending 
on Us, k, T, and /iq such that 

\\he{t)-K^k{t)\\i<Ce^+\ eE (0,eo),tG [0,r']. (5.1) 

(ii) Fix n G N. There is a (3 = l3{s,n) G N such that for any solution ho £ Vp to (jl.3p 
there are constants C and depending onUs,n,T, and ho such that 

||/i,(t)-/i,,„_i(t)|U <C7e", (0,eo),t G [0,r']. (5.2) 

Proof, (i) Let Eq be small enough to ensure that /ie,fc(t) G Us, e £ [0,eo)) t G [0,T]. 
We introduce the differences 

dit) := heit) - he,k{t) and 5{t) := K{e, he{t)) - K{e, he,k{t)). 
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We obviously have 

S{t) = K'{e,The{t) - (1 - T)he,kmd{t)] dr, t £ [0,r'], (5.3) 



and since 

,(l + e2/j/2)3/2 J ' 

we obtain that there exist positive constants ci_2 = ci^2(^s) such that 



K'{e,h)[d\ 



ci\\d\\^<\\5\l.2<C2\\d\l, ae[l,2]. (5.4) 

(Here and in the sequel, we will omit the argument t if no confusion is likely.) In the same 
spirit, for h,h GUg, we introduce the bilinear form B{e, h, h) : i?^(S) x F^(S) M by 

Observe that there are positive constants ci^2 = ci,2(^s) such that 

ci\\d\\l<B{e,h~h){d,d) < C2\\d\\l (5.5) 

as d{t) has zero average over S. 

From (j5.3p and ()3.28p we find, via integration by parts, 

-(9t«i| 5)l2(s) = B{e,he,h^^k){d,dtd) 

= '^(^{dt{B{e,h„h,,k){d,d)) 

- dhB{e,he,hs^k){d, d)dthe - dhB{e,he,he^k)id,d)dths^k^ 

> ^dt{B{e,h„he,k){d,d)) - C\\d\\l (5.6) 
Furthermore, from (|3.3|) and (|4.4|) we have 

dtd{t) =F{e, K{t))Me, Km - F{e, h,,k{t)){K{e, /ie,fc(t))} + R{t) 

=F{e, K{t)){5{t)} + R{t) + R{t), (5.7) 

where 



max||i?(t)||, <Ce^'+i (5.8) 



and 



R:= [ F'{e,Th,^k + {l-r)he)[d]{K{e,h,^k)}dT 
Jo 

By Lemma 13.61 and (|5.4p . 

Plli/2 < C||d||3/2||K(e,/i,,fc)||,„+3/2 < C\\6\\_y2, s G (0,eo), t G [0,T']. (5.9) 
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Multiplying ()5.7p by —6 and applying ()3.19p . (15. 81) . (15.90 . and an interpolation inequality we 
get 

2 _^n\\x\\2 Ce''+M\6\\-i 



-(5td|5)L2(§)<-c||5||f/2_, + C||5||ii/2 + 



< -cMl + (c||<^||g + C\\6f_,) + C (e2'=+2 + ||5||2_^ 

Together with (j5.6p . this shows that 

^B{e, h„ h,^k){d, d) < C(e2^+2 ^ ^^^^ /i„ d)) 
at 

for all e G (0,eo)) ^ £ [0,T']. Taking into consideration that d(0) = 0, we find by Gronwall's 
inequality that 

which proves (|5.ip . 

(ii) Set k := n + bs — 1 and /3 := S2{k). Let eo be small enough to ensure that he^kit) £ i^s, 
ee[0,eo),te[0,T]. 

Instead of ()5.2p we are going to prove the equivalent estimate 

<C7e", EG (0,eo),t G [0,r']. (5.10) 

Let 6s-i := k' {e, he)[d^''~^'>]. Then, in analogy to (IKIHll . 

- (9*^(^-1) I > - C\\d\\l (5.11) 

where iJ(e, /ig) := B{e, hg, hi;). 

On the other hand, differentiating the relation 

dtd = T{e, he) - -F(e, he,k) + R 

{s — 1) times with respect to x we get (cf. Lemma [3^ 

dtd^-^) = -F{e, K)[6s-i] + Ps{e, K, K,k) + R^''^^ ■ 

Recalling ()3.19p . (I3.34p . (15. 4p . and ()5.8|) we obtain from this by Young's inequality 

-{dtd^'-^) I 5s-i)m^^ < -c||5,_i||2/2^, + Ce~^d\U,/2\\6,.i\\y2 + Ce'=+ip,_i||_i 

< -c\\6s-i\\l + Ce-'\\d\\l^,/, + Ce''^' + 

Consequently, by ()5.1ip . ()5.ip . and an interpolation inequality, 

9<i3(e,/ie)((i(^-^\d(^-^))<-c|^f||2+i + Ce-5M||f;i'/^||d||i/^ + Ce^^ 

< C(e-i°^||d||? + e2fc+2) < Ce2fe+2-10s < ^^2n_ 

Integrating over t and using (j5.5p . we obtain (|5.10p . □ 
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Proof of Theorem \1.2[ Choose k and /3 = f3{s,n) as in the proof of Proposition 15.11 (ii), 
let h* := ho{0) and let a and M be such that /io([0, T]) C Ug. By compactness, ^ := 
dist(9Z//s, /io([0,T])) > 0. Let eo be small enough to ensure that he ki[0,T]) C Us, 

distidUs,he,n-ii[0,T])) > /i/2, e G [0,eo), (5.12) 

and Ceq < /i/4, where C is the constant from ()5.2p . 
Let e G (0,eo) and let 

K G c([o,r,),F^-i(S)) nci([o,r,),/7/^-4(S)) 

be a maximal solution to (]3.3p with /ie(0) = h* . In view of Proposition 15.11 (ii), it remains 
to show that > T. Assume < T. The blowup result in Theorem 11.11 (iii) implies that 
there is a T' e (0,r) such that he{[0,T']) C but dist{dUs,h{T')) < fi/A. In view of (fO]) 
and (j5.12p . this is a contradiction to our choice of eo- D 
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